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ABSTRACT 

The study of the angular and spatial structure of the X-ray sky has been under investigation since 
the times of the Einstein X-ray Observatory. This topic has fascinated more than two generations of 
scientists and slowly unveiled an unexpected scenario regarding the consequences of the angular and 
spatial distribution of X-ray sources. It was first established from the clustering of sources making 
the CXB that the source spatial distribution resembles that of optical QSO. It then it became evident 
that the distribution of X-ray AGN in the Universe was strongly reflecting that of Dark Matter. In 
particular one of the key result is that X-ray AGN are hosted by Dark Matter Halos of mass similar 
to that of galaxy groups. This result, together with model predictions, has lead to the hypothesis 
that galaxy mergers may constitute the main AGN triggering mechanism. However detailed analysis 
of observational data, acquired with modern telescopes, and the use of the new Halo Occupation 
formalism has revealed that the triggering of an AGN could also be attributed to phenomena like 
tidal disruption or disk instability, and to galaxy evolution. This paper reviews results from 1988 to 
2011 in the field of X-ray selected AGN clustering. 

Subject headings: Surveys - Galaxies: active - X-rays: general - Cosmology: Large-scale structure of 
Universe - Dark Matter 



1. INTRODUCTION 

After about 50 years from the opening of the X-ray 
window on the Universe with the discovery of Sco-Xl 
and the Cosmic X-ray background (CXB, Giacconi et 
al. 1962), our knowledge of high energy processes in the 
Universe has dramatically improved. One of the lead- 
ing mechanism for the production of X-ray in the Uni- 
verse is accretion onto compact objects. For this rea- 
son the study of astrophysical X-ray sources is a power- 
ful tool for studying matter under the effects of extreme 
gravity. As the efficiency of converting matter into en- 
ergy in accretion processes is proportional to the "com- 
pactness" of the object, (i.e. oc M/R), it is clear that 
the strongest sources powered by accretion are Super- 
Massive Black Holes (SMBH). It also became a corner- 
stone of astrophysics that every galaxy with a bulge-like 
component hosts a SMBH at its centre and that the BH 
mass and the bulge velocity dispersion are strictly re- 
lated (|Magorrian et al.l |1998|) . It is also believed that 
black holes reach those high masses via one or more 
phases of intense accretion activity and therefore shin- 
ing as Active Galactic Nuclei (AGN). It is believed that 
an AGN is basically shine mostly from the power emit- 
ted by a thin, viscous, accretion d isk orbiting the central 
SMBH S hakura fc Sunvaevl (|1976l ). Such a disk produces 
a high amount X-rays both from its hot inner regions (as 
far as the soft X-ray emission is concerned) and from a 
non thermal source which is supposed to be the primary 
source of X-rays (both soft and hard) . 

Since its discovery, the nature of the CXB has 

1 INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, 
40127 Bologna, Italy 

2 University of Maryland, Baltimore County, 1000 Hilltop Cir- 
cle Baltimore, MD 21250, USA 

3 Max-Planck-Institut fur Plasmaphysik and Excellence Clus- 
ter Universe, Boltzmannstrasse 2, D-85748 Garching, Germany 

4 Max-Planck-Institute fur Extraterrestrische Physik, Giessen- 
bachstrasse 1, D-85748 Garching, Germany 



been strongly debated, but soon the community con- 
verged into interpreting most of the CXB as the in- 
tegrated emission of AGN across the cosmic time. 
While the discrete nature of the CXB h as been pro- 
posed (Bergamini Londrillo fc Setti I fl967h and rapidly 
unvei led by experiments li ke Einstein ([Giacconi et all 
Il979h and ROSAT (see e.g. lHasinger et al.lll993[ ). little 
cosmological information has been obtained from sam- 
ples of AGN because of the scarce number of detected 
sources in the X-ray band. Structure formation models 
and numerical simulations have shown that structures 
in the Universe have undergone a hierarchical growth 
starting from the denser peaks in the primordial gaussian 
matter distribution. The Large Scale Structures (LSS) 
of the Universe are gravitationally dominated by Dark 
Matter (DM) and we can consider it as the responsible 
and one of the main driver of the Cosmological struc- 
tures evolution. Dark matter is believed to clump in 
large scale halos (DMH INavarro. Frenk fc White] 119971) 
which are populated by galaxies. Thus galaxies can be 
considered as tracers of the DM distribution in the Uni- 
verse and the study of their spatial clustering led us to 
a most comprehensive view of the LSS. On the other 
hand AGN/Quasar, as phase of the galactic evolution, is 
a quite rare phenomenon in the Universe as their space 
density of these objects is about 1/100-1/1000 lower than 
that of galaxies. This means that AGN/Quasar survey 
require large field of view and/or deep exposure to pro- 
vide statistically significant samples. 
The study of their clustering and its evolution is a pow- 
erful tool to understand, from a statistical point of view, 
what kind of environment is more likely to host AGN. 
This is not just an academic question but, this is strictly 
related to the mechanism of AGN activation. We know 
that one of the candidate me chanism for triggering an 
AGN is galaxy merger (see e.g.lHopkins et al.ll2007ll2008t 
IHopkins fc Henquistl [2001 ISilverman et al.l I201l[) . The 
probability of such an event is definitely dependent on 
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the environment inhabited by the host galaxy. Even if 
the mean distance between galaxies is relatively small, in 
high density (mass) environments, they have a high ve- 
locity dispersion and therefore, the likelihood of a major 
merger is very low. On the contrary, in the field the like- 
lihood of galaxy mergers is low because of the large aver- 
age distance between galaxies. The most favorable place 
to detect a merger is therefore a modera tely low density 
(mas s) environment like a group (see e.g [Mcintos h et all 
[2009]) . 

In fa ct, merger-driven models (see e.g. IHopkins et al.l 
|2QQ7t) accurately predict the observed large-scale clus- 
tering of quasars as a function of redshift up to z— 
4. The clustering is precisely that predicted for small 
group halos in which major mergers of gas-rich galaxies 
should proceed most efficiently. Thus it is well estab- 
lished empirically and with theoretical predictions that 
quasar clustering traces a characteristic host halo mass 
~4x 10 12 ft. _1 Af©, supporting the scenario in which ma- 
jor mergers dominate the bright quasar populations. 

In addition other phenomena like secular processes 
ma y become dominant a t lower l uminosities as suggeste d 
bv iMilosavlievic et al.1 (12006ft : IHopkins etafl (j2006h : 
IHopkins fc HenquistJ ([20091 ) . Low- luminosity AGN could 
be triggered in more common nonmerger events, like 
stochastic encounters of the black holes and molecular 
clouds, tidal disruption or disk instability. This leads to 
the expectation of a characteristic transition to merger- 
induced fueling around the traditional quasar-Seyfert lu- 
minosity division (growth of BH masses above/below 
~ 1O 7 M ). However the triggering mechanism of the 
SMBH growth must be compliant with Mbh-p relation, 



that links the growth of the SMBH with growt h of the 
bulge of the host galaxy (IMagorrian et al.lll998l ). 
As shown in IHopkins et al.l ([20081 ) . the predicted large- 
scale bias of quasars triggered by secular processes is, at 
all redshifts, lower than the bias estimated for quasars fu- 
eled by major mergers. This implies that low-luminosity 
Seyfert galaxies live in DMHs that never rich the char- 
acteristic mass associated with small group scales. 
On the other hand, the majority of the results on the 
clustering of X-ray selected AGN, suggest a picture where 
moderate-luminosity AGN live in massive DMHs (12. 5< 
logMuuH [h~ 1 M ( 7 ) ] <13.5) up to z ~ 2, i.e. X-ray se- 
lected AGN samples appear to cluster more strongly than 
bright quasars. The reason for this is not completely 
clear but several studies argued that these large bias and 
DMH masses could suggest a different AGN triggering 
mechanism respect to bright quasars characterized by 
galaxy merger-induced fueling. 

This paper reviews results of clustering of X-ray se- 
lected AGN from the first Einstein to the most recent 
Chandra and XMM-Newton surveys. We give a detailed 
description of the methods used in this kind of analy- 
sis from simple power-law to halo models. In addition 
we discuss the results of X-ray AGN clustering in the 
framework of AGN evolution and triggering. We adopt a 
ACDM Cosmology with O A =0.7, O n=0.3, H n =100 h~* 
km/s/Mpc with h=0.7 and <t 8 =0.8 ((Larson et al.ll201lL 
WMAP-7). 



PREVIOUS MEASURES OF X-RAY CLUSTERING 
AMPLITUDE 



As far as the X-ray source clustering results are con- 
cerned, the development of the field has always be driven 
by with the performance of the telescopes. In particular 
while first results studied the angular distribution of the 
unresolved CXB under the assumption that Quasars were 
its main contributors, recent Chandra and XMM-Newton 
surveys sample clustering of AGN with a precision com- 
parable to that achievable with redshift galaxy surveys. 
In the following section we will use the following con- 
vention for reporting results of clustering analysis in 
the case of power-law representation of the auto(cross)- 
correlation function: if the clustering is measured in the 
angular space, we will use: 



w{0) = 0/9 Q 



1-7 



(1) 



where 0q is the angular correlation length. If the mea- 
surements has been performed in the real (redshift) space 
this becomes: 

£(r) = r I V ~ 7 (£(s) = s/s ~ 7 ,in z - space), (2) 

where 7 is the 3D c orrelation slope and rn or So are the 
correlation lengths. iBarcons fe F abian (1988) measured 
with Einstein a clustering signal of the CXB on scales 

< 5 corresponding to an angular correlation length 

#0 ~ 4 . They have shown the importance of studying 
the angular structure of the CXB by pointing out that 
a large fraction of the CXB could have been attributed 
to sources with a redshift distribution similar to optical 
QSOs. In additon, the first prediction was not consis- 
tent with the hypothesis that the CXB was also partly 
produced by a diffuse hot Intergalactic Medium (IGM) 
component. It was also proposed that these sources 
were actually clustered on comoving scales of the order 
of ~1Q h' 1 Mpc. 

iCarrera fe Barco"nsl (Il992fl. iGeorgantopoulos et al.l 
(|1993D and iSoltan fe Hasingerl (|1994t) observed that 



the CXB was highly isotropic on scales of the order of 
2°-25°. The first attempt of measuring t he clusterin, 
of X-r ay selected AGN was performed by iBovle fc M< 
(1993), that measured a barely significant signal by 
using a sample of 183 EMSS sources, mostly local AGN 
(z<0.2). These evidences have brought the attention 
to the study of the clustering of the CXB down to the 
arcminute scale. The first significant upward turn for 
the measurement of AGN clustering in the X-ray band 
has been brought to light by ROSAT. By using a set 
of ROSAT-PSPC Pointing on an area of -40 deg 2 , 
iVikhlinin fc Formanl ([19951 ) measured, for the first time 
an angular correlation signal of faint (ROSAT) X-ray 
sources on scales <10 / . By usin g the Limber equation 
(see Appendix and Peebles 1980) they have de- projected 
their angular correlation function into a real space cor- 
relation function and found that, under the assumption 
that the the redshift distribution of the sources was the 
same as that of optical QSOs, the spatial correlation 
length was in the range 6-10 h _1 Mpc. With such a 
result, they confirmed the hypothesis that the CXB was 
mostly produced by sources with a redshift distribution 
comparable to that of optically selected QSO, though 
wi th almost double sou r ce den sity. By using the r esults 
of IVikhlinin fc Formanl (fl995h and lAkvlas et all (2000. 
who obtained similar results), IBarcons et al.l ( 20011) has 
shown for the first time that X-ray selected AGN are 
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highly biased tracers of the underlying LSS at z<l by 
showing a redshift evolving bias factor as large as b~2. 

However, it is worth to consider that the deprojection 
of the angular correlation function into a 3D correlation 
relies on several assumptions, like the model dependent 
expected redshift distribution, which may lead to a bi- 
ased estimate of the real space clustering. It is how- 
ever worth noticing that angular correlation can be very 
useful to provide a first overview in the early phase of 
surveys, when optical identifications are not available, 
especially sampling new part of the parameter space of 
sources, like i.e. new unexplored luminosity /flux limits 
and therefore source classes. Detailed physical models 
are however much better investigated by more sophisti- 
cated techniques as shown in the following parts. 

The first firm detection of 3-D spatial clustering of 
X-ray selected AGN has been claimed by iMullis et al.1 
(|2004h by using data of the ROSAT-NEP survey. They 
detected on an area of ~81 deg 2 a 3a significant signal 
in the redshift space auto-correlation function of soft X- 
ray selected sources at < z >^ 0.22. They have shown 
that at that redshift AGN cluster with a typical correla- 
tion length ro=7.4±1.9 h -1 Mpc. Their results suggest 
that the population of AGN in such a sample is con- 
sistent with an unbiased population with respect to the 
underlying matter . Their result suggested that at that 
redshift AGN were hosted in DMHs of mass of the order 
of 10 13 h' 1 M Q . 

With the development of Chandra and XMM-Newton 
surveys and thanks to the high source surface densities 
(i.e. >400-1000 deg -2 ) our capabilities in tracing the 
LSS has dramatically increased. One of the first evi- 
dences that AGN are highly correl ated with the under - 
lyin g LSS has been pointed out bv lCappi et al. I (|2001h 
and ICappelluti et~afl (|2005l) and references therein, who 
showed that around massive high-z galaxy clusters the 
source surface density of Chandra point sources is signif- 
icantly, up to two times, higher than that of th e back- 
ground. More recently, iKoulouridis fc Plionisl (|2010f ) 
showed that although the X-ray source surface density of 
AGN around galaxy clusters is larger than in the back- 
ground, the amplitude of their overdensities is about 4 
times lower than that of galaxies in the same fields. This 
has been interpreted as a clear indication of an environ- 
mental influence on the AGN a ctivity. Silverman et al. 
(2011) in the COSMOS field and lKoss et all (f2010h in the 
Swift-B AT all-sky survey have shown that the AGN frac- 
tion in galaxy pairs is higher relative to isolated galaxies 
of similar stellar mass providing an additional evidence 
of the influence of the environment on AGN activity. 

Chandra and XMM-Newton performed several blanck 
sky extragalactic surveys, and most of them dedicated 
part of their efforts in the study of the LSS tra c ed by 
AGN to unveil their co-evolution. [Basilakos et al.l (|2004l 
12005( 1 by using data of the XMM-Newton 2dF-survey 
have measured an unexpected high correlation length 

both in the angular (9q ~10 ) and, by projection, in 
the real space (r ~16 h -1 Mpc). Such an high correla- 
tion length has been detected in this field only, thus one 
can explain such a measureme nt as a statistical fluc tua- 
tion. With the same technique, iGandhi et al.l (|2006l ) ob- 
tained a marginal 2-3er detection of angular clustering in 
the XMM-LSS survey and obtained 6> =6.3(42)±3(±[ 3 ) 



in the 0.5-2 (2-10) keV bands and a slope 7 -2.2. 
iPuccetti et al.l (|2006l) measured the clustering of X-ray 
sources in the XMM-Newton ELAIS-S1 survey in the soft 
and hard energy bands with a sample of 448 sources. 

They obtained 9 = 5.2±3.8 a" and 6» =12.8±7.8 4" in 
the two bands respectively. These measurements have 
been deprojected with the Limber's inversion in the real 
space and obtained ro=9.8-12.8 h _1 Mpc and ro=13.4- 
17.9 h _1 Mpc in the t wo bands, re s pectiv ely. 

In the Chandra era, iGilli et all (|2005l ) measured the 
real space auto-correlation function of point sources in 
the CDFS-CDFN. They have measured in the CDFS 
r =8.6±1.2 h -1 Mpc at z=0.73, while in the CDFN they 
obtained ro=4.2±0.4 h -1 Mpc. The discrepancy of these 
measurements has been explained with variance intro- 
duced by the relatively small field of view and the con- 
sequent random sampling o f LSSs in the field. In the 
CLASXS survey lYang et all (|2006l ) obtained a measure- 
ment of the clustering at z=0.94 with ro=8.1^2' 2 n_1 
Mpc which proposes that AGN are hosted by DMH of 
mass of 10 121 h -1 M Q (see next Section). In addi- 
tion they proposed that AGN clustering evolves with 
luminosity and they found that the bias factor evolves 
with the redshift. Such a behavior is similar to that 
found in optically selected quasars. The XMM-Ne wton 
(lHasinger et al.l 120071: ICappelluti et "all 120071 120091) and 



Chandra (jElvis et all 120071 : IPuccetti et al.ll2009f ) survey 
of the COSMOS field have provided a leap forward to 
the field of X-ray AGN clustering by surveying a 2 deg 2 
field of view. The key of th e success of this project is a 
redshift survey zCOSMOS (jLillv et al.ll2007t) performed 
simultaneously with the X-ray survey, together with ob- 
servations in more than 30 energy bands from radio to 
X-ray, that allowed to measure either the spectroscopic 
or the photometric redshift of every source. In the X-ray 
band, the survey covers 2 deg 2 with XMM-Newton with 
a depth of ~60 ks with the addition of a central 0.9 deg 2 
observed by Chandra with ~150 ks exposure. The first 
sample of ^1500 X-ray sources (ICappelluti et al.ll2007l ) 
has been used bv lMivaii et al.l (|2007t ) to determine their 
angular correlation function, without knowing their dis- 
tance, and just assuming a theoretical redshift distribu- 
tion for the purpose of Limber's deprojection. Significant 
positive signals have been detected in the 0.5-2 band, in 

the angular range of 0.5 -24 , while the positive signals 
were at the ~ 2a and 3a levels in the 2-4.5 and 4.5-10 keV 
bands, respectively. With power-law fits to the ACFs 
without the integral constraint term, they have found 

correlation lengths of 6q=1.9 ±0.3 , 0.8 , and 6 

// 

±2 for the three bands, respectively, for a fixed slope 
7=1.8. The inferred comoving correlation lengths were 
r =9.8±0.7, 5.8±{$, and 12±2 h _1 Mpc at the effective 
redshifts of z=l.l, 0.9, and 0.6, respectively. Comparing 
the inferred rms fluctuations of the spatial distribution of 
AGNs as. agn (see Appendix) with those of the underly- 
ing dark matter, the bias parameters of the X-ray source 
clustering at these effective redshifts were found in the 
range b = 1.5-4. Such a result lead to the conclusion that 
the typical mass of the DMH hosting an AGN is of the 
order Mdmh ~10 13 M ^ h" 1 Similar results have been 
found by lEbrero et all p009l ) using the angular correla- 
tion function of 30000 X-ray sources in the AXIS survey. 
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In the XMM-LSS survey lElviv et all (|2011l) measured 
the clustering of ^5000 AGN and computed via Lim- 
ber's deprojection the obtained ro= 7.2±0.8 Mpc/h and 
r = 10.1±0.8 Mpc/h and 7 ~2 in the 0.5-2 keV and 2-10 
keV energy bands, re spectively. In the XMM-COSMOS 
field iGilli et al.1 <|2009f ) measured the clustering of 562 X- 
ray selected and spectroscopically confirmed AGN. They 
have obtained that the correlation length of these source, 
r =8.6±0.5 h- 1 Mpc and slope of 7 =1.88±0.07. They 
also found that if source in redshift spikes are removed 
the correlation length decreases to about 5-6 h _1 Mpc 
. Even if not conclusively, they also showed that nar- 
row line AGN and broad line AGN cluster in the same 
way, indicating that both class of sources share the same 
environment, an argument in favor of the unified AGN 
model which predicts that obscuration, and therefore the 
Type-I/Type II dichotomy is simply a geometrical prob- 
lem. However it is worth noticing that such a procedure 
may artificially reduce the clustering signal and the ef- 
fects of such a cut in the sample, may lead to an unreli- 
able estimate of the cl ustering signal. 

Even if the results of IGilli et al.l <|2009f ) provide a quite 
complete overvi ew of the environment s of the AGN in the 
COSMOS field. lAllevato et al.l (poTTh analyzed the same 
field by using the halo model formalism (see Section 
Their results show that AGN selected in the X-ray band 
are more biased than the more luminous optically se- 
lected QSO. This observation significantly deviates from 
the prediction of model s of merger driven AGN activ- 
ity ([Hopkins et al.l 120061 : iBonoli et al.l I2009D . indicating 
that other mechanisms like disk/bar instability of tidal 
disruptions may trigger an AGN. They also found that 
Type 1 AGN are more biased than Type 2 AGN up to 
redshift of ~ 1.5. 

In the Bootes field [Hickox et al. (2009) explored the con- 
nection between different classes of AGN and the evolu- 
tion of their host galaxies, by deriving host galaxy prop- 
erties, clustering, and Eddington ratios of AGN selected 
in the radio, X-ray, and infrared (IR) wavebands from 
the wide-field (9 deg 2 ) Bootes survey. They noticed that 
radio and X-ray AGNs reside in relatively large DMHs 
(M DM H ~3x10 13 and 10 13 M h" 1 , respectively) and 
are found in galaxies with red and green colors. In con- 
trast, IR AGN are in less luminous galaxies, have higher 
Eddington ratios, and reside in halos with Mdmh <10 12 
M q h" 1 . 

On the same line, iCoil et all (|2009f ) measured the clus- 
tering of non-quasar X-ray active galactic nuclei at z = 
0.7-1.4 in the AEGIS field. Using the cross-correlation 
of Chandra-selected AGN with 5000 DEEP2 galaxies 
they measured a correlation length of ro = 5.95±0.90 
h _1 Mpc and slope 7= 1.66±0.22. They also con- 
cluded that X-ray AGN have a similar clustering am- 
plitude as red, quiescent and "green" transition galaxies 
at z ~1 and are significantly more clustered than blue, 
star- forming galaxies. In addition they proposed a "se- 
quence" of X-ray AGN clustering, where its strength is 
primarily determined by the host galaxy color; AGN in 
red host galaxies are significantly more clustered than 
AGN in blue host galaxies, with a relative bias that is 
similar to that of red to blue DEEP2 galaxies. They 
did not observe any dependence of clustering on opti- 
cal brightness, X-ray luminosity, or hardness ratio. In 



addition they obtained evidence that galaxies hosting 
X-ray AGN are more likely to reside in groups and 
more massive DMHs than galaxies of t he same color 



and luminosity without a n X-ray AGN. lAllevato et"al 
(poll . ICoil et all (pOOl ): iMountrichas fc Georgakakis 



(|2011f ) concluded that DEEP2 X-ray AGN at z - 1 are 
more clustered than optically selected quasars (with a 
2.6cr significance) and therefore may reside in more mas- 
sive DMHs. In an evolutionary picture their results are 
consistent with galaxies undergoing a quasar phase while 
in the blue cloud before settling on the red sequence with 
a lower-luminosity X-ray AGN, i f they are similar ob jects 
at different evo lutionary stage s dHickox et al.ll2009D . At 
lower re dshift, iKrumpe et all (|2010T ) confirmed the re- 
sults of ICoil et all (|2009h . Various recent works have 
been presented indications and/or evidences, of varying 
significance, regarding a correlation between the X-ray 
Luminosity and the AGN cluster i ng amplitude, based ei- 
ther on th e spat ial (lYang et al.l 120061: IGilli et all 120091: 



Coil et all 120091: ICappcllut fet al.l 120101: IKrumpe et "ail 



20101 120111 ) ~ot the angular (jPhonis et all 120081) corre 



lation function. 
Note that luminosity dependent clustering is one of the 
key features of merger triggered AGN activity and is one 
of the prime motivations for AGN clustering analyses. 
Low Lx AGN have been found to cluster in a similar 
way as blue star forming galaxies while high L^ AGN 
cluster like re d passive galaxies . Such a result has been 
confirmed by ICappelluti et al.l ()2010| ) using the Swift- 
BAT all sky survey at z~0. They detected both a Lx 
dependence of AGN clustering amplitude and a larger 
clustering of Type I AGN than that of Type II AGN. 
IKrumpe et al.l (2010, 2011) confirm the weak dependence 
of the clustering strength on AGN X-ray luminosity at a 
2cr level for z < 0.5. 

Table [D] summarizes all the discussed results on 
the clustering of AGN in X-ray surveys with bias 
factors converted to a common cosmology (17 a = 
0.7,O„ =0.3, cr s =0.8) in the EMSS , [Boyle fc Md (fl99l ; 
RASS. lVikhlinin fc Formanl (11995I) : lAkv las et al.l 12 000): 



ROSAT-NEP, iMullis et al.l (12001: A X IS. lEbrero et all 
ll2"009h: E LAIS-S1. IPuccetti et all (120061) : CDFS 
IGilli et all (12001 : CDFN. IGilli et all (120051) : lYang et all 
XMM- 2dF. IBasilakos et all (12001: XM M-LSS. 

(120061) : 
Ml): 



(2006); 



CLAS X S. lYangetal 
(200§; lAllevato et a! 



I Gandhi et alj 
COSMOS. IGilli et al _ 

Swift- BAT, ICappelluti et al.l (120101): A EGIS, lCoiLet_al| 
(120091): AG ES IHl ^ kox et aLI I2009h: ROSAT-SDSS, 
IKrumpe et al.1 (|2010f ). while fig. |3] shows the redshift 
evolution of the correlation length ro as estimated in pre- 
vious works, according to the legend. 

2.1. Techniques of investigation 

The continuously increasing volume and quality of 
data, allowed a parallel improvement of the techniques 
of investigation. The fi rst surveys of Einstein (see e.g. 
iBarcons fc Fabiar]|1988fl . used the autocorrelation func- 
tion of the unresolved CXB and linked it to the clustering 
properties of the clustering of X-ray source that produced 
it. 

Modern surveys have mostly estimated correlation 
function with estimators that use of random samples and 
real data pairs and then estimating physical clustering 
properties by fitting the correlation function functions 
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with simple power-law models in the form of eq. [5] A 
detailed description of the method to estimate correla- 
tion functions is given in the appendix. Considering its 
power, here we give a detailed description of halo mod- 
eling which is by far the most reliable formalism to de- 
scribe clustering of AGN/Galaxies and to determine the 
environment of a specific DMH tracer. 

3. HALO MODEL 

In the hierarchical model of cosmological structure for- 
mation, galaxies, group of galaxies, clusters an so on are 
built from the same initial perturbation in the underlying 
dark matter density field. Regions of dark matter denser- 
than-average collapse to form halos in which structures 
form. Galaxies and AGN, as well as, groups and clusters 
are believed to populate the collapsed DMHs. 

The theoretical understanding of galaxy cluster- 
ing has been greatly enhanced through the frame- 
work o f the halo model (iKauffmann. Nusser fc Steinmet j 



of AGN: 



19971: [Peacqck fe Smitbl 120001: ICoorav fe Shethl I2002t 
Tinke r et all 120051 : iZheng et all 120*051 1 One can fill 
DMHs with objects based on a statistical halo occupa- 
tion distribution (HOD), allowing one to model the clus- 
tering of galaxies within halos (and thus at non-linear 
scales) while providing a self-consistent determination of 
the bias at linear scales. Similarly the problem of dis- 
cussing the abundance and spatial distribution of AGN 
can be reduced to studying how they populate their host 
halos. 

The HOD analysis recasts AGN clustering measurements 
into a form that is more physically informative and con- 
ducive for testing galaxy/AGN formation theories. 

Thus, one can use measurements of AGN two-point 
correlation functions to constrain the HOD of different 
sets of AGN and gain information on the nature of DMH 
in which they live. In fact, the power of the HOD mod- 
eling is the capability to transform data on AGN pair 
counts at small-scales into a physical relation between 
AGN and DMH at the level of individual halos. 

The key ingredient needed to describe the clustering 
properties of AGN is their halo occupation distribution 
function P/v(Af^), which gives the probability of finding 
N AGN within a single halo as a function of the halo 
mass, Mh- In the most general case, P^iMh) is entirely 
specified by all its moments which, in principle, could 
be observationally determined by studying AGN clus- 
tering at any order. Regrettably AGN are so rare that 
their two-point function is already poorly determined, so 
that it is not possible to accurately measure higher-order 
statistics. One overcomes this problem by assuming a 
predefined functional form for the lowest-order moments 
of Pjv(Af^), defining the halo occupation number N{Mh) 
which is the mean value of the halo occupation distri- 
bution N(M h ) =<N > (M h ) = J2 N NP N (M h ). It is 
convenient to describe N{Mh) in terms of a few param- 
eters whose values will then be constrained by the data. 

An accurate description of matter clustering on the ba- 
sis of the halo approach requires three major ingredients: 
the halo mass function n{Mh) (the number of DMHs 
per unit mass and volume), the mass-dependent bias- 
ing factor b(Mh) and the density profile of halos. These 
terms, along with a parametrization of N(Mh), allow us 
to calculate some useful quantities; the number density 



riAGN = / n{M h )N{M h )dM h 



the large-scale bias: 

= / b h (M h )N(M h )n(M h )dM h 
J N(M h )n(M h )dM h 

and the average mass of the host dark halo: 



M 



f M h N{M h )n{M h )dM h 
J N(M h )n(M h )dM h 



(3) 



(4) 



(5) 



The number density and clustering properties of 
the DMHs can be easily computed, at any red- 
shift, by means of a set of analytical tools which 
have been tested and calibrated against numerical 
simulations (iMo fc White! 119961: ISheth fe Tormenl 119991; 
' Sheth et all hoOlt iTinker et al l 120051: jBasilakos et al l 
20081: iTinker et al l[2oTot iPillepich et alH201d : iMa et al.l 
20111 ). Popular choices for both n(Mh) and b(Mh) are 
the analytical spherical collapse ( She th fc Tormeniri999[ ) 
or an ellipsoidal collapse model (ISheth et al.ll2001l see $1] 
for more details) . A detailed description of HOD math- 
ematical formalism is given in Appendix B. 

3.1. Occupation Number 

In the past ten years, a very successful framework 
for modelling the nonlinear clustering properties 
of galaxies has been developed and a number of 
halo models have been presented in the literature. 
These have been successfully used to describe the 
abundance an d clustering pr o pertie s of g alaxies 
at both low dPeacock fc Smith] boOOt iSeliak 
Scoccimarro et al. I 120011: iBerlind fc Weinberg 



2000; 



2002 



Marinoni fc Hudson! I2002t iMagliocchetti fc Porcianil 



2003; 



van den Bosch, Yang fc MoT 



Yang. Mo fc van den Bosch! 12002 
Tinker et al.1 120051: iPhleps et al.1 



Zch avi et al.l 1200 
20061: 



Zheng 



et al.1 



2002 



2004; 



2009 ) and high dBullock. Wcchsler fc Somerville 
Moustakas & So mervilld 120021: lHamana et al 
IZhengl 120041 ; IZheng et al.l 12007! ) redshifts, as well as 
whether these galaxies occupy the centers of the 
DMH or are sat ellite galaxies (|Kravtsov et al.l 120041 : 
IZheng et al.1l2005T ). 

Partially due to the low number density of AGN, 
there have been few results in the literature interpreting 
AGN correlation function using HOD modelling, where 
the small-scale clustering measurements are essential. 
iPorciani. Magliocchetti fc Norberg I ((20041 ) studied the 
clustering of 2QZ QSO with the halo model to infer the 
mean number of optically selected quasars which are har- 
boured by a virialized h alo of given mass and t he ch arac- 
teristic quasar lifetime. iPadmanabhan et "all (|2009f) dis- 
cussed qualitative HOD constraints on their L RG-optical 
QSO Cross- Correlation Function (CCF) and lShen et al.l 
(2010) modelled with the HOD the observed two-point 
correlation function of 15 binary quasar at z > 2.9. 

The standard halo approach used for quasars and 
galaxies is based on the idea that the elements of HOD 
can be effectively decomposed into two components, sep- 
arately describing the properties of central and satellite 
galaxies within the DMH. A simple parametric form used 
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to describe the galaxy HOD is to model the mean occu- 
pation number for central galaxies as a step function, 



i.e.. 

(N c 



(N c , 



1 for halos with mass M > M„ 



and 



for M < M„ 



while the distribution of 



satellite objects can be well approximated by a Pois- 
son distribution with the mean following a power law, 
(N sat ) = (M/Mi) a . Previously derived HOD of galax- 
ies show a values ~ 1 — 1.2 which implies a number of 
satellite galaxies approximately proportional to Mh- 

The clustering properties of X-ray selected AGN have 
been modelled with the HOD in two previous works 
for sources in the Bootes field IStarikova et al.l (120101) 
and in the ROSAT A ll-Sky Survey iMivaii et al.1 (|2011[ ). 
IStarikova et al.l (|2010l ) used the the projections of the 
two-points correlation function both on the sky plane and 
in the line of sight to show that Chandra/Bootes AGN 
are located at the center of DM halos with M > M m i n = 
4 x 10 12 h -1 Mq, assuming a halo occupation described 
by a step function (zero AGN per halo/subhalo below 
M m i n and one above it). They also showed that Chan- 
dra/Bootes AGNs are located at the centers of DMHs, 
limiting the fraction of AGN in non-central galaxies to be 
< 0.09 at the 95% CL. The central locations of the AGN 
host galaxies are expected in the merger trigger model 
because mergers of equally-si zed galaxies preferen tially 
o ccur at the centers o f DMH (|Hopkins et al.ll2008D . 

IMivaii et ail (120111 ) modelled the AGN HOD testing 
the effects of having or not AGN in central galaxies by 
using the RASS AGN-LRG cross-correlation. In the first 
scenario they assumed that all the AGN are satellites 
and they visualized the HOD of the LRG as a step func- 
tion with a step at logM h [h~ 1 M Q = 13.5]. While for- 
mally they assumed that all AGN are not in central 
galaxies, the HOD constraints obtained from this as- 
sumption can be applied to satellite and central AGN 
if the AGN activity in central galaxies of high-mass ha- 
los (logM/ l [h~ 1 M Q > 13.5]) is suppressed. In particular, 
they used a truncated power-law satellite HOD, with 
two parameters: the critical DMH mass below which 
the AGN HOD is zero and the slope a of the HOD for 
Mh > M cr . They also investigated a model where the 
central HOD is constant and the satellite HOD has a 
power-law form, both at masses above M m i n . In all the 
cases they rejected a ~ 1, finding a marginal prefer- 
ence for an AGN fraction among satellite galaxies which 
decreases with increasing Mh- They argued that this re- 
sult might be explained by a decrease of the cross-section 
for galaxy merging in the environment of richer groups 
or clusters. In fact previous observations infer that 
the AGN fraction is smaller in cluster s than in groups 
([Arnold et al.ll2009HMartini et alj|2009t iSilverman et al.1 
2009; Ko ulouridis fc Plionisl^oToj ). 

It is important to stress that the small number statis- 
tics has so far limited the accuracy of correlation func- 
tion of X-ray AGN at small-scales, especially through 
the auto-correlation function of the AGN themselves. 
The situation can be improved by measuring the cross- 
correlation function of AGN with a galaxy sample that 
has a much higher space density, with common sky and 
redshift coverage as the AGN redshift surveys. The 
AGN clustering through cross-correlatio n function with 



galaxies is emerg i ng in the last years 
Coil et al. I l2007t ICoil et al.1 120091 Ifflc 



(ILi et al 



Hic kox et al.1 2009 ; 
Padmana bhan et al.ll2009t iMountrichas et all 120091) and 



2006 



Press & Schechter (1974) 
Sheth & Tormen (1999) /' 
Sheth et al. (2001) / 
Tinker et al. (2005) / 



/ 



/ 
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/ 



logM„ [M s h-'] 

Fig. 1. — Halo bias as function of halo masses for a fixed red- 
shift z=l and the corresponding predictions of Pr ess & Shechter 
1974) (long- dash ed line), [Sheth & Tormen (19991, dashed line) 

■ ■ otz ■ 



Sheth et al T (t200ll . solid line) and lTinker et ah! (120051 . dotted line). 

can be used to i mprove our under s tanding of how AGN 
populate DMH (jMivaii et al.ll201lHKrumpe et al.ll2011l) . 

4. BIAS AND DMH MASS 

In the literature, the bi as parameter is often calculated 
with the power-law fits (IMullis et al.l 120041 lYang et al. | 
[20061: [Mivaii et al.l|2007l: ICoil et al.l 120091: iKrumpe et al.l 
l201fi Iffickox et al.l 120111 ) over scales of 0.1-0.3 < r p < 
10-20 ft." 1 Mpc. The power-law models of the ACF are 
usually converted to the rms fluctuation over 8 ft. -1 Mpc 
spheres or are averaged up to the distance of 20 ft -1 Mpc. 
While some authors use only large scales (r p > 1 — 2 ft -1 
Mpc) to ensure that the linear regime is used, others 
include s maller scales to ha ve better statistics. As an 
example, Hick ox et al.l (2009) fitted their data with a bi- 
ased DMH projected correlation function. 
In the HOD analysis the bias factor only comes fr o m the 
2-halo term (r p >l-2 ft -1 Mpc). IMivaii et all (|2011h 
compared the bias of RASS-AGN from the full HOD 
model (Eq. ID4[) with the one estimated using the power- 
law best fits parameters, finding that the bias estimates 
are consistent within la. Moreov er, using Eq . ID1I one 
introduces large statistical errors. lAllevato et al.1 (2011) 
found a similar results in comparing the bias of X-ray 
AGN in COSMOS field from the 2-halo term with Eq. 
ID3I and the one estimated from the power-law best fits 
parameters. In Appendix C, we describe the mathemat- 
ical procedures for the bias parameter calculation com- 
monly used in the lite rature. 

Most of the authors (iHickox et al.ll2009l : IKrumpe et al.1 
12011 ICappelluti et al.l 1201017 " us ed an analytical expres 
sion (as the one described in 



Sheth et all 120011 : iTinker et al.l 120051: iBasilakos et al 
20081) " to assign a characteristic DMH mass to the host 



Sheth fc f orTner] Il99 



1 



ing halos. The large-scale bias is directly related to the 
mass function of halos, so that the mass of a halo dictates 
the halo clustering and the number of such halos. The 
halo mass can be quantified in terms of the peak height 
v = S c /o~(Mh, z), which characterizes the amplitude of 
density fluctuations from which a halo of mass Mh form 
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at a given redshift. In general one assumes 8 C = 1.686 
and a(Mh,z) is the linear overdensity variance in spheres 
enclosing a mean mass Mh- The traditional choice of 
the mass function and th en of the bias has been that of 
IPress fc Schechterl (Il97l : 



b PS = l 



v 2 -\ 



(G) 



A commonly-used p rescription was derived by 
ISheth fc Tormenl (fl99S 



lST 



2 P /d c 



1 + {av 2 )P 



(7) 



where a = 0.707 and p = 0.3 or the ellipsoidal collapse 
formula of ISheth et all (l200ll) : 



<SMT 



av 



\fab{c 



(av 2 ) c 



+6(l-c)(l-c/2)] 



(8) 



where a = 0.707, b = 0.5, c = 0.6 or the recali- 
brated parameters a — 0.707, b — 0.35, c = 0.8 of 
I Tinker et al.l ([2005). The v parameter can be est imated 
following the Appendix of Ivan den Boschl (|200l . Fig. 
[T] shows the bias as function of t he halo mass Mh, 
at z== l, following the pre d iction s of IPress fc Schechterl 
(fl97l: ISheth fc Tormenl (fl999h : ISheth et al.l (|2001[k 
iTinker et al.l (I2005D. 

lAllevato et al.l ( 20111) argued that this approach reveals 
an incongruity due to the fact that the AGN bias used in 
the formulas above, is the average bias of a given AGN 
sample at a given redshift. In fact, following this ap- 
proach one can not take into account that the the average 
bias is sensitive to the entirety of the mass distribution; 
different mass distributions with different average masses 
can give rise to the same average bias. 

On the contrary by using the halo model, the average 
bias and the average mass of the sample, Eq. ID4l and Eq. 
Improperly account for the shape of the mass distribution: 
the average bias depends on the halo number density and 
on the AGN HOD, integrated over the mass range of the 
particular AGN sample. They introduced a new method 
that uses the 2-halo term in estimating the AGN bias 
factor assuming an AGN HOD described by a (5-function. 
Following this approach they properly took into account 
for the sample variance and the growth of the structures 
over time associated with the use of large redshift interval 
of the AGN sample. 

On the o t her h and, iMivaji et all (|2011| ) and 
iKrumpe et al.1 (|2011h applied the HOD modeling tech- 
nique to the RASS AGN-LRG CCF in order to move 
beyond determining the typical DMH mass based on the 
clustering signal strength and instead constrain the full 
distribution of AGN as a function of DMH mass. Along 
with a parametrization of N(Mh) they estimated the 
large-scale bias and the typical mass of hosting DM ha- 
los using Eq. ID4I and Eq. [5] This method improves 
the clustering analysis because it properly uses the non- 
linear growth of matter in the 1-halo term through the 
formation and growth of DMHs. These results are signif- 
icant improvements with respect to the standard method 



of fitting the signal with a phenomenological power law 
or using the 2-halo term (see Appendix C). 

4.1. X-ray selected AGN bias, bias evolution and mass 
of the hosting halos 

The majority of the X-ray surveys agree with a pic- 
ture where X-ray AGN are typically hosted in DM 
halos with mass of the order of 12. 5< logMpMH 
[h^M p] <13.5 at low (z<0.4) an d high (z~l) red- 
shift (iGilli et al.l [2005k lYang et all 1200a iGilli et al 



20091: Iffickox etalJl2009f ICoil et al.ll2009b IKrumpe et al 
2010: ICappelluti et al.l 120101: iStarikova et al.l 120101: 
Mivaii et al.ll2011h IKrumpe et alj|2011l) . 



At high redshift, IGilli et al.l (|2005f) measured the clus- 
tering of X-ray AGN with z=0-4 in both the ~ 0.1 deg 2 

for 240 sources in the north- 



1.87 



+0.14 
-0.16 



2.64^q3q for 124 sources in the southern 



CDFs, finding b 
ern field and b = 
field. At z ~ 1. lYang et al.l (|2006l ) measured the cluster- 
ing of 233 spectroscopic sources in the 0.4 deg 2 Chandra 
CLASXS area and of 252 spectroscopic sources from the 
CDFN, both at z=0.1-3. They found b = 3.58±?;|| for 
the CLASXS AGN and b = 1.77±g*fg for the CDFN field. 
IGilli et al.1 <f2009h studied 538 XMM-COSMOS AGN 
with 0.1 < z < 3 and they found a bias factor b — 
3.081QJ4 at I ~ 1. Using the Millennium simula- 
tions they suggested that XMM-COSMOS AG N reside in 
DMH with mass M DMH > 2.5 x lO 12 /*" 1 M . ICoil et al.l 
(2009) measured the clustering of X-ray AGN at z=0.7- 
1.4 in the AEGIS field a nd th ey estimated b = 1.85±g;||. 
Following Zhe ng et al.l (|2007l ) they infer from the bias 
factor that at z = 0.94 the minimum DM halo mass of the 
X-ra y AGN is > lO^M^fe" 1 . These r esults combined 
withlMountric has fc Ge orgakakis (2011) show that mod- 
erate luminosity X-ray selected AGN live in DMHs with 
masses Mh ~ lO 13 ^ 1 ./^ at all red shifts since z ~ 1 . At 
lower redshift iHickox et all (|2009h analysed 362 AGES 
X-ray AGN at < z >=0.51. The bias factor equal to 
b = 1.40 ±0.16 indicates that X-ray AGN inhabit DM 
halos of typical mass ~ 10 13 M^h -1 . 

In the local Universe lCappelluti et al.1 (|2010l ) estimated 
for - 200 Swift-BAT AGN a bias equal to b = 1.21±g;g£ 
which corresponds to logM£>M=13.15 lo.i3 n_1 -^0- 
lAllevato et al.l ()201lD estimated an average mass of 
the XMM-COSMOS AGN hosting halos equal to 
logM [h- 1 Mpc] = 13.10 ±0.06 at z ~ 1.2. They also 
measured the bias of Type 1 and Type 2 AGN, find- 
ing that the latter reside in less ma ssive halos than 
Type 1 AGN. Only two other works (|Cappelluti et al.l 
120101 : IKrumpe et al.ll20Tfj|) analysed the clustering prop- 
erties of X-ray se l ected Type 1 AGN and Type 2 AGN. 
ICappelluti et al.l (|2010l ) estimated the typical DM halo 
mass hosting type 1 and type 2 Swift-BAT AGN at 
z 0. They measured that these two different sam- 
ples are characterized by halos with mass equal to 
logMDMih^Mo] - 13.94±g;|f and - 12.92±g:||, respec- 
tively. However the lack of small separation pair of Type 
I AGN in the local Universe may have produced system- 
atic deviations which were not accounted in their fits. In 
IKrumpe et all (l2010h the bias factor of BL RASS AGN 
at z = 0.27 are consistent with BL AGN residing in halos 
with mass logM DM [h,- 1 M & ] = 12.58±g;|g. 
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Fig. 2. — Bias factor (Left Panel) and mass of AGN hosting halos (Right Panel) as a function of redshift for X-ray selected AGN (black 
data points), X-ray selected Type 1 AGN (blue data points) and X-r ay selected Type 2 AGN (red data p oints) as esti mated in different 
survey s (COS MOS, IGilli et ail 120091) ; lATTevato et al.l i20W) -, CDFN, Gilli ct al. (2005); Yan g et al.l 12009); Swift-BAT , Cappellut i et al.l 
120101): CDFS.IGilli et al.l 1)2003): AEGIS.ICoil et al.l (|2009l l: AGES. Iffickox et al.l potisT ): ROSAT-NEP. IMullis et all pOOl ); ROSAT-SDSS, 
Krumpc ct al. (2010); CLASXS, Yang ct al. (2006)). The dashed lines show the expected b(z) of DMHs with different masses according to 
the legend, based on Sheth et al. (2001). The grey points show results from quasar - quasar correlation measurements using spectroscopic 

samp les from SDSS IIRoss et a.l.|[2509t ISTieh et al.ll2509T> ■ 2QZ (ICroom et al.ll20"05l ; [Porciani fc Norberg|[2006T) and 2SLAQ (Ida Angela et al.l 
120081 ). All the previous studies infer the picture that X-ray selected AGN which are moderate luminosity AGN compared to bright quasars, 
inhabit more massive DMHs than optically selected quasars in the range z = 0.5 — 2.25. 

Using the HOD model, IStarikova et al.l ([20101 ) suggested tions with SDSS galaxies and derive the bias by applying 



that X-ray Chandra/Bootes AGN are located at the cen- 
ter of DM ha l os wit h M > M min = 4 x 10 12 /i _1 M Q while 
iMivaji et all ([20111 ) estimated for RASS AGN at z=0.25 
b = 1.32 ± 0.08 and a typical mass of the host halos of 
13.09 ±0.08. 

The redshift evolution of the clust ering of X- r ay se - 
lected AGN has been first studied bv lYang et all (|2006f ) 
in the CLAXS+CDFN fields. They measured an increase 
of the bias factor with redshift, from b = 0.95 ± 0.15 at 
z=0.45 to b — 3.03 ±0.83 at z=2.07, corresponding to an 
a verage halo mas s of ~ 12.11 h _1 M©. 

lAllevato et al.l ([201 ll ) studied the redshift evolution 
of the bias for a sample of XMM-COSMOS AGN at 
z < 2. They found a bias evolution with time from 
b(z = 0.92) = 1.80 ± 0.19 to b(z = 1.94) = 2.63 ± 
0.21 with a DM halo mass consistent with being con- 
stant at logM [ft _1 M ] ~ 13.1 at all redshifts z < 2. 
They also found evidence of a redshift evolution of the 
bias factor of XMM-COSMOS Type 1 AGN and Type 
2. The bias evolves with redshift at constant average 
halo mass logM [h- 1 M Q ] ~ 13.3 for Type 1 AGN and 
logMoih^Mol ~ 13 for Type 2 AGN at z < 2 2 5 and 
z < 1.5, respectively. In particular I Allevato et all ([20111 ) 
argued that X-ray selected Type 1 AGN reside in more 
massive DMHs compared to X-ray selected Type 2 AGN 
at all redshifts at ~ 2.5a level, suggesting that the AGN 
activity is a mass triggered phenomenon and that differ- 
ent AGN classes are associated with the DM halo mass, 
ir respective of redshift z . 

IKrumpe et al.l ([201 ID measured the clustering ampli- 
tudes of both X-ray RASS and optically-selected SDSS 
broad-line AGNs, as well as for X-ray selected narrow- 
line RASS/SDSS AGNs through cross-correlation func- 



thc HOD model directly to the CCFs. They estimated 
typical DMH masses of broad-line AGNs in the range 
logiMh/l^MQ]) = 12.4 - 13.4, consistent with the halo 
mass range of typical non-AGN galaxies at low redshifts 
and they found no significant difference between the clus- 
tering of X-ray selected narrow-line AGNs and broad- line 
AGNs up to z ~ 0.5. 

Fig. [5] shows the bias parameter (Left Panel) and 
the mass of the AGN hosting halos (Right Panel) as a 
function of redshift for X-ray selected AGN (black data 
points), X-ray selected Type 1 AGN (blue data points) 
and X-ray selected Type 2 AGN (red data points) as esti- 
mated for different surveys (see the legend) . The dashed 
lines show the expected b(z) of typical DM halo masses 
Mdmh based on Sheth et al. (2001). The masses are 
given in logM dmh in units of h^ 1 M Q . 

There have been several studies of the bias evolu- 
tion of optical quasar with the redshift as shown in 
fig- 121 (grey data points), based on large survey sam- 
ples such as 2QZ, 2SLAQ and SDSS (ICroom et al.ll20F 
IPorciani fc NorbereJ [2001 IShen et al.l 120091: IRoss et al 
2009t lda Angela et al.ll2008[ ). These previous studies in- 
fer the picture that X-ray selected AGN which are moder- 
ate luminosity AGN compared to bright quasars, inhabit 
more massive DMHs than optically selected quasars in 
the range z = 0.5 — 2.25. 

Recently, IKrumpe et al.l ([20111 ) verified that the clus- 
tering properties between X-ray and optically- selected 
AGN samples are not significantly different in three 
redshift bins below z — 0.5 (the differences are 1.5<r, 
O.lcr and 2. Oct). The reason for the fact that X-ray se- 
lected AGN samples appear to cluster mor e strongly than 
optically- selected AGNs is still unclear. Allevato et al. 



AGN Clustering 

(|2011h ; iMountrichas fc Georgakakisl (|201lD suggested 
that the difference in the bias and then in the host 
DMH masses is due to the different fueling mode of those 
sources from that of the X-r ay selected moderate lurni- 
nosity AGN. On the contrary, iKrumpe' "elaTl (pnTTh sug- 
gested that some of the X-ray clustering studies signifi- 
cantly underestimate their systematic uncertainties and 
then it may turn out that these measurements are consis- 
tent with optical AGN clustering measurements. More 
high-z AGN clustering measurements based on larger 
samples are needed to gain a clearer picture. 

4.2. AGN Life Time 

One of the most important tests for studying the evo- 
lution models of AGN is understanding their lifetime. It 
is widely accepted that AGN is phase of the galaxy life 
necessary to explain the coevolution of the bulge and the 
black hole. After a triggering event, of which we do not 
know the nature, yet, the central black holes begins its 
accretion phase and it is believed that undergoes several 
regimes of Eddingto n rates and bolometric luminosity. 
iMartini fc Weinberg] (|2001) proposed a method to de- 
rive the AGN life time by knowing their space density 
and their DMH host mass. 

By knowing the AGN and DMH halo space density 
at a given luminosity and mass (xiagn, ^dmh), one 
can estimate the duty cycle of the AGN, tagn(z) — 
n AGN (L z) r \ Tffoste r(H(z)) is the Hubble time at 
a given redshift. Actually this method provides only an 
upper limit since it assumes that the life of an halo of a 
given mass is similar to the Hubble time. A more exhaus- 
tive formulation would t agn (z) = n*^(M% TDMH ( z ) > 
where tdmh{z) is the age of a DMH at given redshift. 
Unfortunately this quantity cannot be estimated ana- 
lytically but could be estimate in a statistical way by 
using hidrodynamic simulations. Several results can be 
mentioned for this quantities but their dispersion is very 
large, ther e fore w e report only some example. At z=l, 
IGilli et all (|2009D obtains t hat the typical duty cyc le of 
AGN is <1 Gyr. At z=0, ICappelluti etaLl (pOlOl ) has 
measured a duty cycle in the range 0.2 Gyr- 5 Gyr with 
an expectation value of 0.7 Gyr. Both the measurements 
are fairly larger than the 40 million years determined by 
IMartini fc Weinberg] (|2001[ ) at z=2-3. These differences 
however are not surprising if we assume that the differ- 
ent populations of AGN, grow with a different Eddington 
rate a s function of their t ypical luminosities and/or red- 
shifts (jFabian et al.ll2009ft . 

5. DISCUSSION 

In this paper we reviewed the results in the field of 
X-ray AGN clustering, for energies between 0.1 keV to 
55 keV over a period of more that 20 years. The liter- 
ature has produced an increasingly convincing and con- 
sistent picture of the physical quantities derivable from 
this kind of study. Most of the advancements in the field 
have been achieved with the improvement of survey ca- 
pabilities and instruments sensitivity. The availability of 
simultaneously wide and deep fields, coupled with multi- 
wavelength information, has produced larger and larger 
samples of spectroscopically confirmed sources. This al- 
lowed several teams to refine the techniques needed to 
estimate the two point ACF and the quantities derived 
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form it. In particular we are entering a phase where, 
at least at z< 2, AGN clustering studies won't probably 
provide any new result unless evaluated with the HOD 
formalism. Open questions as what is the AGN occu- 
pation number and the evolution of HOD define a new 
barrier which is necessary to break in order to under- 
stand the history of X-ray emission from accretion onto 
AGN. In this respect, samples of X-ray selected AGN al- 
ways need a spectroscopical follow-up to provide a solid 
base to compute clustering in the real space rather than 
in the angular space. 

Summarizing, the current picture is that X-ray selected 
AGN are highly biased objects with respect to the un- 
derlined matter distribution. Such an evidence is clearer 
when measuring the redshift dependence of AGN bias. 
At every redshift from z=0 to z=2, AGN cluster in way 
similar to DMH of mass of the order of log(M Q /i _1 )=13. 
The spread of such a value is of the order 0.3-0.5 dex 
at ler. This means that the determination of what kind 
of environment is inhabited by AGN, is relatively well 
constrained and identical at every redshift sampled by 
X-ray surveys. This allows us to formulate the hypoth- 
esis that every phase of AGN activity is mass triggered 
phenomenon (i.e. each AGN evolutionary phase is char- 
acterized by a critical halo mass). 

It is believed that major mergers of galaxies is one of 
the dominant mechanisms for fueling quasars at high red- 
shift and bright luminosities, while minor interactions, 
bar instabilities or tidal disruptions are important at 
low re dshift (z < 1) and low luminosities (L < 10 44 er<7 
s" 1 ) dHopkins et all [20M [200l iHasinger et all l200l 
Hopk ins fc Henquistl I2009I K In the local Universe, for 
example, the st udy of the envir onment of Swift BAT 
Seyfert galaxies (|Koss et al. 2010) finds a larger fraction 
of BAT AGNs with disturbed morphologies or in close 
physical pairs (<30 kpc) compared to matched control 
galaxies or optically selected AGNs. The high rate of 
apparent mergers (25%) suggests that AGN activity and 
merging are critically linked for the moderate luminos- 
ity AGN in the BAT sample. Moreover models of ma- 
jor mergers appear to naturally produce many observed 
properties of quasars, as the quasar luminosity density, 
the shape and the evolution of the quasar luminosity 
function and the large- scale quasar clustering as a func- 
tion of L and z (e g iHqpkins et al l [2001 iShenl [200l 
Shan kaF et alj 120091 120101: iShankarl l2010t iBonoli et al.l 
2009; Trci ster et al. 1 120101) . Quasar clustering at all red- 
shift is consistent with halo masses similar to group 
scales, where the combination of low velocity dispersion 
and moderate galaxy space densit y yields to the high- 
est probability of a cl ose encounter ([Hopkins et al.ll2008l : 
iMcIntosh et al.ir2009f ). Moreover recent detections of an 
hx dependent clustering play in favor of major mergers 
being the dominant AGN triggering mechanism. 

On the other hand it has became clear that many AGN 
are not fueled by major mergers and only a small fraction 
of AGN are associated with mo rphologi cally disturbed 
galaxi es. iGeorgakakis et all ()20071 ) and ISilverman et all 
( 2009) found that AGN span a broad range of environ- 
ments, from the field to massive groups and thus ma- 
jor mergers of galaxies, possibly relevant for the more 
luminous quasar phenomenon, may not be the primary 
mechanism for fueling these moderate luminosity AGN. 
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iGeorgakakis et all (|2009l) suggest that bar instabilities 
and minor interactions are more efficient in producing lu- 
minous AGN at z < 1 and not only Seyfert galax i es and 
low-luminosity AGN as the Hopkins fc Henquistl (2009) 
model predicts. ICisternas et al.l (|2010t) analysed a sam- 
ple of X-ray selected AGN host galaxies and a matched 
control sample of inactive galaxies in the COSMOS field. 
They found that mergers and interactions involving AGN 
hosts are not dominant and occur no more frequently 
than for inactive galaxies. Over 55% of the studied AGN 



sample that is characterized by Lbol ~ 10 45 erg s _1 and 
by mass of the host galaxies M* > 1O 1O M0 are hosted by 
disk-dominated galaxies, suggesting that secular fuelling 
mechanisms can be highly efficient. 

Moreover several w o rks on the AGN h ost galaxies 
(IDunlop et al.l 120031: IGrogin et all 120051; IPierce et al.| 
20071: IGabor et al.l 120091; IReichard et al. 2009JMT et al.l 

2009) show that the morphologies of the AGN host galax- 



ies do not present a preference for merging systems. 

At high redshift (z ~ 2) recent fin dings of 
iSchleeel et al~l (|200lD and lRosario et al.l (|2011|) , who ex- 
amined a smaller sample of AGN in the ERS-II region 
of the GOODS-South field, inferred that late-type mor- 
phologies are prevalent among the AGN hosts. The role 
that major galaxy mergers play in triggering AGN activ- 
ity at 1. 5 < z < 2.5 was also st udied in the CDF-S. At 
z= 1.5-3 Schawinsk i et al.l ()201lD showed that for X-ray- 
selected AGN in the Chandra Deep Field South and with 
typical luminosities of 10 42 erg s _1 <Lx <10 44 erg s _1 
the majority ( 80%) of the host galaxies of these AGNs 
have low Srsic indices indicative of disk-dominated light 
profiles, suggesting that secular processes govern a signif- 
icant fraction of the cosmic growth of black holes. That 
is, many black holes in the present-day universe grew 
much of their mass in disk-dominated galaxies and not 
in early- type galaxies or majo r mergers. 

Later. iKocevski et al.l (|2011h found that X-ray selected 
AGN at z ~ 2 do not exhibit a significant excess of dis- 
torted morphologies while a large fraction reside in late- 
type galaxies. They also suggested that these late- type 
galaxies are fueled by the stochastic accretion of cold gas, 
possibly triggered by a disk instability or minor interac- 
tion; 

lAllevato et al.l (|201lD argued that for moderate lumi- 
nosity X-ray AGN secular processes such as tidal disrup- 
tions or disk instabilities might play a much larger role 
than major mergers up to z ~ 2.2. 

It becomes important to study the clustering properties 
of AGN at high redshift when we assume the peak of 
the merger-driven accretion. Moreover given the com- 
plexity of AGN triggering, a proper selection of AGN 
samples, according to the luminosity or the mass of the 
host galaxies can help to test a particular model boost- 
ing the fraction of AGN host galaxies associated with 
morphologically disturbed galaxies. 

From the evolutionary point of view the evidence of 
a bias segregation of optically and X-ray selected AGN 
might be a sufficient proof to claim that the two phenom- 
ena are sensitive to different environments and therefore 
likely driven by different triggering mechanisms. A more 
comprehensive picture will be available when the clus- 
tering of different phases of AGN activity will be studied 
a nd compared. 

Hick ox et al.l (|2009l ) interpreted their clustering results 
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Fig. 3. — Redshift evolution of the correlati on length r p as es - 
timated in diff e rent X -ray survey s (COSMO S . IGilli et al.1 120091'); 
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12003'! ; Swift-B AT. ICappelluti et al. (2010); CDF S. IGilli et al.l 
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in terms of a general picture for AGN and galaxy evolu- 
tion which is reproduced in Fig. 01 The picture consists 
of an evolutionary sequence that occurs at different red- 
shifts for halos with different masses. In this scenario, 
luminous AGN accretion occurs preferentially (through 
a merger or some secular process) when a host DMH 
reaches a critical Mdmh between 10 12 and 10 13 M Q 
h (this phase is indicated by the solid ovals). Once 
a large halo reaches this critical mass, it becomes visi- 
ble as a ULIRG or SMG (owing to a burst of dusty star 
formation) or (perhaps subsequently) as a luminous, un- 
obscured quasar. The ULIRG/quasar phase is associated 
with rapid growth of the SMBH and formation of a stel- 
lar spheroid, and is followed by the rapid quenching of 
star formation in the galaxy. Subsequently, the young 
stellar population in the galaxy ages (producing " green" 
host galaxy), and the galaxy experiences declining nu- 
clear accretion that may be associated with an X-ray 
AGN. Eventually the aging of the young stars leaves a 
"red" and "dead" early- type galaxy, which experiences 
intermittent "radio- mode" AGN outbursts that heat the 
surrounding medium. For "medium" initial DMHs, the 
quasar phase and formation of the spheroid occurs later 
than for the systems with high halo mass, so that at 
z ~ 0.5 we may observe the green X-ray AGN phase. 
Even smaller halos never reach the threshold mass for 
quasar triggering; these still contain star-forming disk 
galaxies at z < 0.8, and we observe some of them as op- 
tical or IR-selected Seyfert galaxies. The dashed box in- 
dicates the AGN types (in their characteristic DMH) that 
would be observable in the redshift range 0.25 < z < 0.8. 
Further steps in the field will require the study of clus- 
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Fig. 4. — Schematic for a simple picture of AGN and host galaxy evolution, taken from Hickox ct al. (2009), and motivated by the AGN 
host galaxy and clustering results presented in that study. 



tering of AGN from z=3 to z=6-7. This will likely lead 
to the determination of the mass of early DM spheroids 
who hosted primordial black holes seeds. However this is 
a very challenging task since it requires a very deep and 
wide surve y with an almost co mplet e optical follow-up. 

BOSS (lEisenstein et all I20H and BigBOSS 
(jSchlegel et al. 1120011 ) wiiTdetect high redshift AGNs at 
z ~ 2.2, which will improve AGN clustering measure- 
ments at higher redshifts. The only approved mission 
that at the moment wi ll allow to study the z=3-5 X-ray 
Universe is eROSITA (jPredehl et all 120071 launch Dec. 
2013) for which an estimate of the completeness of the 
typical follow up is still unavailable. Additional l y, the 
Large Synoptic Survey Telescope (jlvezic et al.l 120081 
LSST) is expected to identify ^2 million AGNs in op- 
tical bands. eROSITA and LSST have the potential to 
significantly improve AGN clustering measurements at 
low and high redshifts, though only if there are dedicated 
large spectroscopic follow-up programs. Another strong 
contribution will come form either Nustar that will likely 
provide a better view of AGN clustering without the 
selection biases introduced by photoelectric absorption. 
Athena the proposed ESA new generation telescope that 
will mount a wide field imager on a very large collecting 
area telescope, will provide a further view on the deep 



X-ray sky and likely push our knowledge of the high-z 
X-ray Universe. 

In addition to better model the evolution of SMBH 
environments a fundamental point to start is to establish 
the nature of BH seeds at z=10. Such a determination 
will likely come with the new generation of telescope like 
JWST and ESO-ELT. 

NC thanks the INAF-Fellowship program for sup- 
port. NC thanks the Delia Riccia and Blanccflor- 
Lodovisi-Boncompagni foundation for partial support. 
VA is supported by the DFG cluster of excellence Origin 
and Structure of the Universe (www.universe-clustcr.de). 
NC, VA and AF thank the referees, Ryan Hickox and 
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the paper. 
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APPENDIX 

DERIVING THE TWO POINTS AUTO-CORRELATION FUNCTION 

The two-point auto-correlation function (£(r), ACF) describes the excess probability over random of finding a 
pair with an object in the volume dVi and another in the volume (LV2, separated by a distance r so that dP = 
n 2 [l + £,(r)]dVidV2, where n is the mean space density. A known effect when measuring pairs separations is that 
the peculiar velocities combined with the Hubble flow may cause a biased estimate of the distance when using the 
spectroscopic redshift. To avoid this effect it is usually computed the projected ACF (| Davis fc PeebleslllQSl : w(r p ) = 
2 Jo™"* t;(r p ,Tr)dn. Where r p is the distance component perpendicular to the line of sight and ir parallel to the line of 
sight (Fisher et al. 1994). It can be demonstrated that, if the ACF is expressed as £(r) = (r/ro) -7 , then 

w(r p ) = A( 7 )ry p -^ (Al) 

where A( 7 ) = r(l/2)r[( 7 - l)/2]/T( 7 /2) (Peebles 1980). 

The ACF is mostly estimated by using the minimum variance estimator described by Landy & Szalay (1993): 



. DD-2DR + RR , A . 

e(rp ' ] = RR (A2) 



where DD, DR and R R a re the normalized number of data-data, data-random, and random-random source pairs, 
respectively. Equation IA2I indicates that an accurate estimate of the distribution function of the random samples is 
crucial in order to obtain a relia ble estimate of £(r p , 7r). Note that other estimator have been proposed in the literature, 
but the lLandv fc Szalay (1993) one has been shown to provide the smallest statistical variance. Such a formalism can 
be easily adopted when computing the angular or the redshift space correlation function, with the only difference that 
the evaluation is made on a single dimension. Several observational biases must be taken into account when generating 
a random sample of objects in a X-ray flux limited survey. In particular, in order to reproduce the selection function of 
the survey, one has to carefully reproduce the space and flux distributions of the sources, since the sensitivity in X-ray 
surveys is not homogeneous on the detector and therefore on the sky. This points out the necessity of create a random 
sample which i nclud es as many selection effects as possible since the estimate of £(r) (or w(0)) is strongly dependent 
on RR (see eq. IA2[) . Moreover in several case optical follow-up of the X-ray source is not 100% complete, therefore 
one must carefully reproduce the mask effect. What is usually done is that to create random samples in 3D, sources 
are placed at the same angular position of the real sources and redshift are randomly drawn from a smoothed redshift 
distribution of the real sources. If instead the spectral completeness is close to 100% then the right procedures is to 
occupy the survey volume with random sources drawn from a L-z dependent luminosity function and accept check if 
they would be observable using a sensitivity map. An important choice for obtaining a reliable estimate of w(r p ), is 
to set n m ax in the calculation of the integral above. One should avoid values of Ti m ax too large since they would add 
noise to the estimate of w(r p ). If, instead, TT max is too small one could not recover all the signal. Uncertainties in the 
ACF are usually evaluated with a bootstrap resampling techniq ue but it is wor t h not ing that in the literature, several 
methods are adopted for errors estimates in two-point statistics iNorberg et a l. (2009, for a detailed description). It is 
known that Poisson estimators generally underestimate the variance because they do consider that points in ACF are 
not statistically independent. Jackknife resampling method, where one divides the survey area in many sub fields and 
iteratively re-computes correlation functions by excluding one sub-field at a time, generally gives a good estimates of 
errors. But it requires that sufficient number of almost statistically i ndependent sub- fields, this is not the case for 
most of X-ray surveys where the source statistics is moderately low. iCoil et all (|2009D estimated the error bars on 
the two-point correlation function including both Poisson and cosmic variance errors estimated, using DEEP2 mock 
catalogs derived from the Millenium Run simulations. 

LIMBER'S DEPROJECTION 

The 2D Angular Correlation Function (ACF) is a projection of the real-space 3D ACF of the sources along the 
line of sight. In the following discussions and thereafter, r is in comoving coordinates. The relation between the 2D 
(angular) ACF and the 3D ACF is expressed by the Limber equation (e.g., Peebles 1980). Under the assumption that 
the scale length of the clustering is much smaller than the distance to the object, this reduces to 
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where dA(z) is the angular distance, N is the total number of sources, and dN/dz is the redshift distribution (per z) 
of the sources. The redshift evolution of the 3D correlation function is customarily expressed by 
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where e=-3 and e—"f- 3 correspond to the case where the correlation length is constant in physical and comoving 
coordinates, respectively. In these notations, the zero-redshift 3-D correlation length ro can be related to the angular 
correlation length 9q by 
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]r(i/2) 

where r(z) is the look-back time. We also define the comoving correlation length 

r o(zeff) = r (l + z- ff )- 3 -^, 



(B4) 



at the effective redshift z e //, which is the median redshift of the contribution to the angular correlation (the integrand 
of the second term) . An essential ingredient of the deprojection process is the redshift distribution of the sources and 
when individual redshifts are not available this is derived from integration of the luminosity function. 

l-HALO AND 2-HALO TERMS IN THE HOD FORMALISM 

In the halo model approach, the two-point correlation function of AGN is the sum of two contributions: the first 
term (1-halo term) is due to the correlation between objects in the same halo and the second term (2-halo term) arises 
because of the correlation between two distinct halos: 

w p (r p ) = w Pt ih(r P ) + w Pt 2h(r P ) (CI) 

Recent articles prefer to express w = (1 + W\h) + w 2 h (|Tinker et all 120051: [Zheng et all 120051: IBlake et al.1120081) . 
instead of w = W\h + Wih, as used in older articles. This is because 1 + £ represents a quantity that is proportional to 
the number of pairs oc [1 + £i/J + [1 + £,2h\- I n this new convention, the projected correlation function w Pt \h represents 
the projection of 1 + rather than 

Similarly, one express the power spectrum of the distribution of the AGN in terms of the 1- and 2-halo term contribu- 
tions: 

P(k) = P lh (k) + P 2h (k) (C2) 

and then the projected correlation function as: 



w P ,ih{r p ) 
w p ,2h{r p ) 



k 

—P lh (k)J (kr p )dk 
k 

—P2h{k)Jo(kr p )dk 
Ztt 



(C3) 
(C4) 



where Jq(x) is the zeroth-order Bessel function of the first kind. 

Several parameterizations exist in literatur e for representing the DMH profile (|Coorav fc Shethl l2002t 
iKnollmann et alJl2008t iStadel et al1l2009[ ) and the iNavarro. Frenk fc White! (|1997j ) (NFW) profile is a popular choice. 
If y(k, Mh) expresses the Fourier transform of the NFW profile of the DMH with mass M/j, normalized such that 
volume integral up to the virial radius is unity, then the one-halo term of the power spectrum can be written as: 

Pih(k) = ^— [ n(M h )N(M h )\y(k,M h )\ 2 dM h (C5) 



n AGN 



Assuming the linear halo bias model (|Mo fc White! fl996l) . the two-halo term of the power spectrum reduces to: 

1 2 



P2h{k)=P m {k) 



1 



nAGN 



n(M h )b(M h )y(k,M h )dM h 



(C6) 



Since the clustering on large scales is do minated by the two-halo ter m, it is fairly insensitive to the assumption of 
AGN distribution inside the hosting halo (|Berlind fc W einberg] 120021 ) . It should be noted that since y ~ 1 on large 
scales (e.g. scales much larger than the virial radius of halos), on such scales the two-halo term can be rewritten as: 



or in terms of projected correlation function: 



P 2h {k)^b 2 P m {k,z) 



v P ,2h(r p ) = b 2 w m , 2h (r p ) 



(C7) 



(C8) 



where b is the bias parameter of the sample and W m> 2h is the DM projected correlation function. For the matter power 
spectrum, P m (k), on e can use the primordial power spectrum with a fixed n s and a transfer fun ction calculated using 
the fitting formula of Eise nstein fc Hul ([T998) or the nonlinear form given bv lSmith et aLI ([2003), Tinker et al. (2005). 
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BIAS PARAMETER CALCULATION 

In t he majority of works on clustering of X-ray AGN (jMullis et al.ll2004t lYang et alJl2006t iGilli et alJl2005t iCoil et al.l 
120091 : IKrumpe et alj|2010t ICappelluti et al. 2010) the standard approaches used to estimate the bias are based on the 
power- law fit parameters of the AGN correlation function. This method assumes that the projected correlation function 
is well fitted by a power-law and the bias factors are derived from the best fit parameters r$ and 7 of the clustering 
signal at large scale. Using the power-law fit one can estimate the AGN bias factor using the power-law best fit 
parameters: 

bpL = cts,,agn{z)/(Tdm{z) (Dl) 

where <Js.agn(z) is the rms fluctuations of the density distribution over the sphere with a comoving radius of 8 Mpc h , 
&dm{z) is the dark matter correlation function evaluated at 8 Mpc fe" 1 , normalize d to a value of <Jdm{ z = 0) = 0.8. 
For a power-law correlation function this value can be calculated by (|PeeblesHl980l ): 

where J 2 (j) = 72/[(3 - 7) (4 - 7) (6 - 7)2^. 

Differently in the halo model approach, the 2-halo term of the projected correlation function, which dominates at 
large scales, can be considered in the regime of linear density fluctuations. In the linear regime, AGN are biased tracers 
of the dark matter distribution and the bias factor is described by: 

b = {w p s h (r p )/w m ,2h(r p )) 1/2 (D3) 

HOD modeling is currently the optimal method to establish the large-scale bias parameter, provided the parametriza- 
tion of N{Mh), by using: 

= Jb h (M h )N(M h )n(M h )dM h 

J N(M h )n(M h )dM h [ ' 

assuming the halo mass function n(Mfe)and the halo bias factor b(Mh). 

In fact, power law fit bias measurements commonly use smaller scales (< 1 — 2h~ 1 Mpc) that are in the 1-halo term 
in order to increase the statistical significance. If power law fits are restricted only to larger scales, the method suffers 
from the problem that the lowest scale, where the linear biasing scheme can still be applied, varies from sample to 
sample and remains ambiguous. 

HOD modeling allows, in principle, the use of the full range of scales since the method first determines the 1 and 2-halo 
terms and then constrains the linear using data down to the smallest r p values that are dominated by the 2-halo term 
for each indivi d ual sa mple. 

IKrumpe et ail (|2011[) estimated the RASS-AGN bias following the power-law (Eq. EU and the HOD (Eq. |D4| 
approach, pointing out that using the first method the err ors on the bias are mu ch larger, but the values are statistically 
consistent which those derived fro m the H OD m odel fits. lAllevato et al.l (|2011l ) found similar results in estimating the 
COSMOS- AGN bias following Eq. iDll and lD3l 

In order to derive a reliable picture of AGN clustering, bias parameters should be inferred from HOD modeling, 
or at least from the comparison of the correlation function with that of the DM only in the linear regime, because 
systematic errors based on power-law bias parameters will be larger than the statistical uncertainties of the clustering 
measurement. 
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